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The electronic and magnetic properties of neutral substitutional transition-metal dopants in dia- 
mond are calculated within density functional theory using the generalized gradient approximation 
to the exchange-correlation potential. Ti and Fe are nonmagnetic, whereas the ground state of 
V, Cr and Mn are magnetic with a spin entirely localized on the magnetic ion. For Co, Ni, and 
Cu, the ground state is magnetic with the spin distributed over the transition-metal ion and the 
nearest-neighbor carbon atoms; furthermore a bound state is found in the gap that originates from 
the hybridization of the Sd-derived level of the dopant and the 2p-derived dangling bonds of the 
nearest-neighbor carbons. A p-d hybridization model is developed in order to describe the origin of 
the magnetic interaction. This model predicts high-spin to low-spin transitions for Ni and Cu under 
compressive strain. 
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I. INTRODUCTION 

Diamond has attracted great interest recently for quan- 
tum information science because of its wide band gap 
and the existence of more than 500 optically-addressible 
centers^, many with long spin coherence times. An 
extensively-studied color center is the NV center, con- 
sisting of a carbon vacancy with a nearest-neighbor (NN) 
substitutional nitrogen atom. In its negatively-charged 
state, the NV~ center carries a spin = 1, which is split 
by the Csv symmetry environment of the spin centeP\ 
its spin state has been proposed as a qubit for quan- 
tum computation. Spin initialization and readout has 
been demonstrated by optical pumping and absorption 
measurement at the zero phonon line (ZPL) at 1.945 
eV. Spin manipulation has been achieved by microwave 
manipulation^ ^, and also through the interaction of the 
N V~ center spin with a nitrogen nuclear spirP or another 
nearby substitutional nitrogen electron spin*". However, 
the NV~ center has a large phonon sideband at room 
temperature that may limit its use for quantum informa- 
tion processing^. 

Other applications, like quantum optics and quan- 
tum cryptography, require stable, bright, and narrow 
linewidth single photon sources (SPS). Potential optical 
centers for SPSs in diamo nd a re the NV~ cente r, t he 
silicon- vacancy (SiV) centei j^n^, Ni-related centers^l^IHIS 
and Cr-related center4^^^. The large phonon sideband 
of the NV~ center spreads the emission over 100 nm re- 
sulting in a Debye- Waller factor of less than 5%, which is 
not ideal for quantum key distributioiP. The SiV center, 
formed by a substitutional Si next to a vacancy, has been 
proposed as a stable SPS with a ZPL at 738 nm, but its 
fast nonradiative decay degrades its fluorescenc^^ and 
ultimately its reliability. 

Transition-metal (TM) dopants have also been studied 
in diamond, however initial studies focused on TM impu- 
rities in natural diamond samples, or in some synthetic 



diamond due to the use of the TM as a solvent cata- 
lyst during the high-pressure-high-temperature (HPHT) 
growth procesJ^^H^. Among the Ni-related centers, the 
NE8 center, formed by a substitutional Ni with four NN 
nitrogen atoms, is responsible for a ZPL at 793 nm with 
a narrow phonon sideband, but its low reproducibility 
impedes its use for applications. Among Cr-related cen- 
ters, a ZPL at 749 nm^^ has been attributed to an in- 
terstitial chromium and two other ZPL at 744 and 764 
nm have been detected but the source has not yet been 
identified^^. TM dopants can also be introduced by ion 
implantatiorP^'I^ or during chemical vapor deposition 
(CVD) -*^^^ , which offer the possibility of introducing the 
dopants substitutionally and without additional atoms 
or as part of a defect complex. Co- and Ni-related cen- 
ters, with different nitrogen, boron and vacancy com- 
plexes, have bee n stu died in various charge states by 
GGA calculationP^J^, finding neutral substitutional Co 
to be stable in pure diamond for a Fermi energy Ep be- 
tween 3.0 and 3.6 eV, and neutral substitutional Ni to be 
stable for between 2.6 and 3.0 eV. 

Transit ion- metal spin centers in diamond offer sev- 
eral potential advantages over other spin centers, includ- 
ing (1) the availability of both strongly- hybridizing (^2- 
symmetry) and weakly-hybridizing (e-symmetry) states 
associated with the d levels^^^^, (2) the high symmetry 
of the substitutional dopant (tetrahedral group, which 
means the crystal field does not split angular-momentum 
1 triplet^^, and (3) the availability of larger spin-orbit 
interactions through the d electrons, which may permit 
faster manipulation of a spin through an electric field 
{e.g. Ref . [26]) . Challenges for these spin centers include 
the large strain introduced when they are placed substi- 
tutionally (although formation energies have previously 
been found to be moderate) , and the potential lower spin 
coherence time due to the larger spin-orbit interaction. 

To clarify these phenomena we present a first- 
principles study of substitutional transition-metal (TMg) 
dopants in diamond. Whereas TM^ dopants may exist 
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in several diflFerent charge states in diamond, we concen- 
trate our study on the neutral substitutional TM dopant 
(TM^) and the chemical trends within the 3d row of the 
periodic table. For each TM^ dopant we determine the 
ground-state configuration, including the charge and spin 
state. We find for Ti through Cu a single p-d hybridiza- 
tion model explains the electronic and magnetic proper- 
ties of the defect. Sc and Zn, which occur at the ends of 
the row, are not described by this model. In Section [TH 
after a description of the method,we present the results 
for 3d transition-metal dopants, in Section [TTI| we describe 
changes under strain of the ground-state spin for Ni and 
Cu and contrast with Cr (which is largely unchanged by 
strain). In Section IV we comment on the implications 
of these results and the strain-induced spin changes for 
quantum information processing. 



II. DOPANT CONFIGURATION 
CALCULATIONS AND p-d HYBRIDIZATION 
MODELS 

A. Method 



and Mn^, which we shall see below become magnetic 
in a similar fashion to Cr^. A different origin is found 
for magnetism in Co^ ,Ni^ and Cu^, which have much 
smaller magnetic energies ~ —100 meV. Sc^ and Zn^ are 
included to complete the 3d row, however the magnetism 
of these ions differs from that of all the other 3d ions. 
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TABLE L TMCes SGGA Results : R is the NN relaxation, 
Em = ^SGGA — Egga the magnetic energy, Mt the total 
magnetization, Mtm (Mc) the on-site magnetization on the 
TM ion (NN C atom), A^tm (A^c) the calculated number of 
electrons on the TM ion (NN C atom). 



The calculations were performed with the scalar 
relativistic version of the full potential local orbital 
FPL09. 00-33 code^*^ using the (spin-polarized) gen- 
eralized gradient approximation ((S)GGA) with the 
parametrization of Perdew, Burke and Ernzerhof^^. The 
convergence of the results with respect to k-space in- 
tegrals was carefully checked. We found that a 8x8x8 
k-point grid were sufficient to obtain convergence of the 
total energy difference between spin-polarized (SGGA) 
and spin unpolarized (GGA) calculations. 

We used a 64-atom TMCes (TM = Sc, Ti, V, Cr, Mn, 
Fe, Co, Ni, Cu, Zn) supercell corresponding to a 2 x 2 x 2 
multiple of the diamond conventional cell with a TM 
substituted for one carbon. We fixed the supercell to 
correspond to the experimental lattice constant^^ of dia- 
mond ao = 3.5668 A, and all the atomic positions within 
the supercell were allowed to relax with a precision of 1 
meV/A. We considered here only T^^-symmetric atomic 
relaxation. As a check, we did a fully-relaxed-atoms cal- 
culation without any symmetry constraints (Ci)for the 
NiCes supercell. The Ci relaxed calculation is about 2 
meV higher in energy than the relaxed calculation, 
which rules out such a symmetry-breaking distortion. 



B. Transition- metal dopant ground-state spin 
configurat ions 

The SGGA results are summarized on Table [H The 
most stable magnetic solution is obtained for Cr^, with 
a magnetic energy Em = —1006 meV, corresponding to 
an = 1 ground state Em lower in energy than the non- 
magnetic GGA solution. Less stable magnetic solutions 
{Em ^ —250 meV) are obtained for the 3d neighbors 




Sc Ti V Cr Mn Fe Co Ni Cu Zn 
FIG. 1. TMCes SGGA total magnetic moment. 



Figure [T] shows the calculated total magnetization of 
TMCes supercells. Ti and Fe are non-magnetic. We 
obtain for Cr and Ni a total magnetization Mt = 2 /^b, 
corresponding to a total spin 5* = 1, although we shall 
see later that these two spin centers differ greatly due 
to the importance of different d orbitals in creating the 
ground-state spin. Zn, also with 5 = 1, differs from both 
as it has a closed d shell. For Sc, V, Mn and Co, we 
obtain a total magnetization Mt = 1 /iB, corresponding 
to a total spin S = ^. The highest total magnetization 
is obtained for Cu with Mt = 3 /iB corresponding to a 
total spin = |. 

For V^, Cr^ and Mn^, the total magnetization of 
Mt = 1.0, 2.0, and 1.0 jub is entirely localized on the 
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TM ion, whereas for Co^, Ni^ and Cu^, the total mag- 
netization of Mt = 1.0, 2.0, and 3.0 fiB is distributed 
between the TM ion and the NN carbon atoms with a 
TM ion magnetization of Mtm = 0.6, 0.9 and 0.8 jub 
and a NN carbon atom magnetization of Mq = 0.1, 0.2 
and 0.4 /iB, corresponding to a total magnetic moment 
0.4, 0.8 and 1.6 /iB distributed over the NN C. The num- 
ber of electrons located on the TM ion compared to its 
atomic number, and the number of electrons on the NN 
carbon atoms allow us to identify the configuration of the 
TM ion and the induced charge on the NN C's. 

From Ti to Fe, the TM occupation number is consistent 
with a TM2+ configuration As^Sd"^ (n=2,3,4,5,6) and the 
TM^+ ion induces two extra electrons distributed over 
the 4 NN carbon dangling bonds with a local charge on 
each NN carbon atoms Qc ^ 0.5. For Co, Ni, Cu, the 
TM occupation number is consistent with a TM+ config- 
uration 45°3(i^+^ (n = 7, 8, 9) and the TM+ ion induces 
one extra electron distributed over the 4 NN carbon dan- 
gling bonds with a local charge on each NN carbon atoms 
Qc ~ 0.2. Sc and Zn differ greatly from all the other 3d 
dopants, and hence will be discussed only at the end of 
Section [m 




Ti V Cr Mn Fe Co Ni Cu 



FIG. 2. Single-particle TM 3d-derived SGGA Kohn-Sham 
(KS) energy levels of TMC63 calculation. TM 3d levels are 
split by the tetrahedral crystal field into e'^ and t2 levels. 
The bonding and antibonding levels (t^ and t^^) are the 
results of the hybridization between the TM t2 and NN carbon 
2p-derived t2 levels. Ey (energy reference) and Ec are the 
SGGA valence band maximum and (direct) conduction band 
minimum of pure diamond. In parenthesis are indicated the 
occupation number of the KS energy levels. 



The TM 3(i-derived Kohn-Sham states within the se- 
ries are presented in Figure (2] The GGA direct and 
indirect band gap for pure diamond are respectively 
^GGA = Ec - Ey = 5.6 and ^^^^^ = 4.3 eV, both of 
which are smaller than the experimental gap Eq^^ = 5.47 
eV due to the well known gap underestimation of the 



GGA functional^. A substitutional TM ion is at the cen- 
ter of a tetrahedron formed by the 4 NN carbon atoms, 
implying a crystal field of symmetry. This crystal 
field will split the single particle 3d manifold into a dou- 
bly degenerate e level and a triply degenerate t2 level. 
For tetrahedral bonding the t2 levels are expected to be 
higher in energy than the e levels. The symmetry of the 
bonds from the four nearest neighbors is compatible with 
t2 symmetr}!^, but not with e symmetry. Thus the 3d- 
derived t2 level hybridizes with the NN carbon 2p-derived 
t2 states to form bonding ts and antibonding Iab hy- 
brid states. Due to the strong hybridization between the 
TM ion and the NN carbon atoms, the 36Z-deriv ed hybrid 
states will form bound states in the gapP^I^^I^. 

For the non-magnetic Ti^ impurity, there are 3 bound 
states in the gap : a totally occupied (6-fold degenerate) 
ts level situated at Ey + 0.6 eV, a totally empty (4- fold 
degenerate) e level at + 3.3 eV and a totally empty 
(6-fold degenerate) Iab level at Ey + 4.8 eV. For V^, 
Cr^, Mn^, the spin up and spin down (3-fold degenerate) 
ts levels are totally occupied and nearly spin degener- 
ate. For V, the Ib forms a bound state at Ey +0.7 eV 
and for Cr^ to Cu^, the ts state forms a crystal field 
resonance state inside the valence band which (with in- 
creasing atomic number) decreases in energy. For V, Cr 
and Mn, the occupation difference of the spin up and spin 
down e bound states is responsible for the total magne- 
tization localized on the TM ion; the Cr e level config- 
uration is (e^)^(e^)^, corresponding to a total magnetic 
moment Mt = 2.0 /j^b {S = 1). The absence of mag- 
netism for Fe is due to a totally occupied spin-degenerate 
e level (4- fold degenerate). For the next three dopants 
in the series (Fig. [T]), Co, Ni and Cu, the e levels remain 
totally occupied and enter the valence band. For these 
dopants the spin up antibonding (3-fold degenerate) 
state starts to become occupied, with 1, 2 or 3 electrons 
respectively, while the spin down antibonding level 
remains totally empty. For Co^, Ni^ and Cu^, the oc- 
cupation difference of the spin up and spin down Iab 
bound states is responsible for the total magnetization 
that is distributed over the TM ion and the NN carbon 
atoms through the hybrid Iab level; for example, the 
Ni tAB level configuration is (^^^)^(^^b)^' correspond- 
ing to a total magnetic moment Mt — 2.0 \ib {S = 1) 
distributed over the Ni ion and the NN carbon atoms. 

A p-d hybridization model^^ helps explain the origin 
of the magnetic interaction and differentiate the behavior 
of those dopants to the left of Fe from those to the right. 
In pure diamond, each carbon shares 4 electrons to form 
covalent bond. If one carbon atom is removed to form a 
vacancy, 4 electrons will occupy the dangling bonds with 
the neutral configuration 2s^2p^ . When a TM ion is in- 
serted at the vacancy site to form a substitutional dopant, 
it will induce an extra charge on the dangling bond dis- 
tributed over the 4 NN carbons. For TM=V,Cr and Mn, 
the magnetization is driven by the e bound states with 
the magnetization entirely localized on the TM ion in 
the TM^+ configuration As^3d^ and two extra electrons 
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shared by the non-magnetic NN carbon atoms. These 
two extra electrons wih occupy the danghng bonds and 
form a defect level of configuration 2s^2p^. The TM 3(i- 
derived t2 levels hybridize with the 2p-derived t2 states 
occupied by 4 electrons. For TM=Co,Ni and Cu, the 
magnetization is driven by the antibonding Iab bound 
states. The TM is in the TM+ configuration 45°3(i^+^ 
which induces one extra electron distributed over the 
NN carbon atoms. This extra electron will occupy the 
dangling bonds and form a defect level of configuration 
2s^2p^ which induces a non-zero magnetic moment on 
the NN carbon atoms. The TM 3(i-derived t2 levels will 
hybridize with these 2p-derived t2 states occupied by 3 
electrons. 

In the next subsections, we will present the electronic 
structure of the different TM^ dopants sorted by their 
spin values. Dopants with the same spin values can have 
dramatically-different origins for the magnetism, depend- 
ing on whether the magnetism is driven by the e states 
or the Iab states. This difference will also lead to a very 
different response of the ground-state spin to strain. 



C. The 5' = centers : Ti° and Fe° 

Fig. [3] shows the total, Ti and NN carbon partial den- 
sity of states (DOS) of the TiCea superceh. The Ti U 
levels are split by the Td crystal field into e and ^2 lev- 
els. Ti t2 levels hybridize with the NN carbon 2p levels 
to form bonding Ib and antibonding Iab hybrid levels. 
The spin-degenerate bonding ts levels are totally occu- 
pied and form a bound state at 0.6 eV above the valence 
band maximum ( VBM) . The spin-degenerate e levels are 
totally empty and form a bound state at 3.3 eV above 
the VBM. The spin-degenerate antibonding Iab levels 
are totally empty and form a bound state situated at 4.7 
eV above the VBM. The number of electrons located on 
the Ti ion is 20.1 which corresponds to a Ti^+ configu- 
ration As^'^d? and spin = 0. Ti^+ induces two extra 
electrons distributed over NN carbon dangling bonds of 
configuration 25^2p^ and spin ^'2 = 0. The Ti (^2)^ level 
hybridizes with the NN carbon 2p^ t2 defect level. The 
spin-degenerate bonding Ib states are totally occupied 
and the e levels are totally empty, which explains the ab- 
sence of magnetism. Fig. [3] explains the observed DOS 
by a p-d hybridization model. The two Ti electrons 
occupy the spin-degenerate Ti ^2 level (violating Hund's 
rule) to form a lower energy closed-shell bonding t b level 
with the four ^2 electrons contributed by the nearest- 
neighbor carbons. 

Fig. [4] shows the total, Fe and NN carbon partial DOS 
of the FeCea supercell. The Fe ?>d levels are split by 
the Td crystal field into e and ^2 levels. Fe ^2 levels hy- 
bridize with the NN carbon 2p levels to form bonding Ib 
and antibonding Iab hybrid levels. The spin-degenerate 
bonding Ib levels are totally occupied and form a crystal 
field resonance state at -1.7 eV below the VBM. The spin- 
degenerate e levels are totally occupied and form a bound 
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FIG. 3. (color online) Top : Total (black solid line), Ti (red 
dashed line) and NN C (blue, dot-dashed line) partial DOS. 
For these and subsequent such plots the upper portion of the 
plot is for the majority spin and the lower portion is for the 
minority spin; here the partial DOS's are the same as the ion 
is nonmagnetic. The bonding state is delocalized across the 
Ti ?>d and NN C 2p states. Bottom : Schematic representation 
of the p-d hybridization model for Ti^ . 



state at 0.9 eV above the VBM. The spin-degenerate an- 
tibonding tAB levels are totally empty and form a bound 
state situated at 3.8 eV above the VBM. The number of 
electrons located on Fe is 24.6, which corresponds to a 
Fe^+ configuration of 4s^3(i^ and spin = 0. Fe^+ in- 
duces two extra electrons distributed over the NN carbon 
dangling bonds, which have a 2s^2p'^ configuration and 
spin S2 = 0. The Ti (^2)^ level hybridizes with the NN 
carbon 2p^ ^2 defect level. The spin-degenerate bonding 
tB and e levels are totally occupied explaining the ab- 
sence of magnetism. Fig. [4] explains the observed DOS 
by a p-d hybridization model, showing it is very similar 
to Ti, but with a full e level instead of the empty e level 
of Ti. 
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FIG. 4. (color online) Top : Total, Fe and NN C partial DOS. 
Bottom : Schematic representation of the p-d hybridization 
model for Fe^. The e states are indicated in red to emphasize 
their localized character; t2 states are indicated in black and 
are delocalized between the Fe and NN C. Symbols are the 
same as in Fig. [3] 



D. The S = h centers : V° Mn° and Co° 
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FIG. 5. (color online) Top : Total, V and NN C partial DOS. 
As before the upper portion of the plot is for the majority 
spin and the lower portion is for the minority spin. Bottom : 
Schematic representation of the p-d hybridization model for 
V°. Symbols are the same as in Fig. 3 



Fig. [5] shows the total, V and NN carbon partial DOS 
of the VCqs supercell. The V 3d levels are split by the 
Td crystal field into e and ^2 levels. V ^2 levels hybridize 
with the NN carbon 2p levels to form bonding and 
antibonding Iab hybrid levels. The bonding ts levels 
are totally occupied and form a nearly spin-degenerate 
bound state at 0.7 eV above the VBM. The spin up and 
spin down e levels are localized on the V site, giving rise 
two bound states situated at 2.5 eV and 3.5 eV above 
the VBM and separated by a Hund exchange splitting 
= 1 eV. The V e level contribution to the total mag- 
netization is about 60 % and the bonding ts states are 
weakly spin-polarized. The antibonding Iab states are 
totally empty and form a nearly spin-degenerate bound 
state at 5.1 eV above the VBM. The number of electrons 
located on the V ion is 21.2, which corresponds to the 
V^+ configuration As^Sd^ and spin Si = ^. V^+ induces 
two extra electrons distributed over the NN carbon dan- 
gling bonds of configuration 25^2p^ and spin S2 = 0. The 



V (^2)^ level hybridizes with the NN carbon 2p^ t2 defect 
level. The bonding Ib states are fully occupied (as for Ti 
and Fe), but in the case of V the e levels are partially oc- 
cupied in the configuration (e^)^(e^)^ corresponding to 
a total magnetic moment Mt = 1.0 fiB entirely local- 
ized on V. Fig. |5] explains the observed DOS with a p-d 
hybridization model. 

Fig. |6] shows the total Mn and NN carbon partial DOS 
of the MnCes supercell. The Mn 3d levels are split by 
the Td crystal field into into e and ^2 levels. Mn ^2 levels 
hybridize with the NN carbon 2p levels to form bonding 
tB and antibonding Iab hybrid levels. The bonding tB 
levels are nearly spin degenerate and form a crystal field 
resonance state located around -2.2 eV below the VBM. 
The spin up and spin down e levels are localized on the 
Mn site, giving rise to two bound states situated at 0.6 
eV and 1.8 eV above the VBM separated by a Hund ex- 
change splitting = 1 eV. The Mn e level contribution 
to the total magnetization is about 70 % and the bonding 



6 



20 



-20 




spin down 



-4 -2 

Energy [eV] 



S 




Cnn DB 



*2t 2p^ 



FIG. 6. (color online) Top : Total, Mn and NN C partial DOS. 
Bottom : Schematic representation of the p-d hybridization 
model for Mn^. Symbols are the same as in Fig.lsl 



tB states are weakly spin-polarized. The antibonding Iab 
levels are weakly spin-polarized and form a bound state 
situated at 5.1 eV above the VBM. The number of elec- 
trons located on the Mn ion is 23.5, which corresponds to 
a 4s^3d^ configuration and spin Si = ^. Mn^+ induces 
two extra electrons distributed over the NN carbon dan- 
gling bonds with 2s^2p^ configuration and spin 6*2 = 0. 
The Mn (^2)^ level hybridizes with the NN carbon 2p^ 
t2 defect level. The bonding Ib states are fully occupied 
and, as in the case of Mn, the e levels are partially oc- 
cupied in the configuration (e^)^(e^)^, corresponding to 
a total magnetic moment Mt = 1.0 /ie entirely localized 
on the Mn. Fig. |6] explains the observed DOS by a p-d 
hybridization model. 

Fig. [7] shows the total, Co and NN carbon partial DOS 
of the CoCes supercell. The Co 3d levels are split by the 
Td crystal field into e and ^2 levels. The Co ^2 levels 
hybridize with the NN carbon 2p levels to form bonding 
tB and antibonding Iab hybrid levels. The bonding tB 
levels are weakly spin-polarized and form a crystal field 
resonance state located around -2.7 eV below the VBM. 
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FIG. 7. (color online) Top : Total, Co and NN C partial DOS. 
The antibonding state tAB has approximately a 50 % Co 3d 
character (50 % NN C 2p character). Bottom : Schematic 
representation of the p-d hybridization model for Co^ . Sym- 
bols are the same as in Fig. |3] 



The spin up and spin down e levels are localized on the 
Co site, giving rise to two peaks at the top of the valence 
band. The number of electrons located on the Co ion 
is 25.9 which corresponds to a Co+ configuration 4:S^3d^ 
and a spin Si = 0. Co+ induces one extra electron dis- 
tributed over NN carbon dangling bonds of configuration 
2s'^2p^ and spin S2 = The Co (^2)^ level hybridizes 
with the NN carbon 2p^ t2 defect level. The 2p^ config- 
uration of the NN carbon dangling bond is responsible 
for the S = ^ spin {Mt = 1.0 /ie) which is distributed 
over the Co {Mqo = 0.6 /ie) and the 4 NN carbon atoms 
(Mc = 0.1 /ie). Fig. [7] explains the observed DOS by a 
p-d hybridization model. 



E. The 5' = 1 centers : Cr° and Ni°. 

Fig. [S] shows the total, Cr and NN carbon partial 
SGGA DOS of the CrCea supercell. The Cr 3d levels are 
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FIG. 8. (color online) Top : Total, Cr and NN C partial DOS. 
Bottom : Schematic representation of the p-d hybridization 
model for Cr^. Symbols are the same as in Fig. 



total magnetic moment Mt = 2.0 /iB entirely localized 
on the Cr ion. Fig. [S] explains the observed DOS by a 
p-d hybridization model. 
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split by the crystal field into e and t2 levels. Cr t2 lev- 
els hybridize with the NN carbon 2p levels to form bond- 
ing ts and antibonding Iab hybrid levels. The bond- 
ing ts levels are nearly spin-degenerate and form crystal 
field resonance states around -1.5 eV below the VBM. 
The spin up and spin down e levels are localized on the 
Cr site, giving rise to two bound states approximately 
0.8 eV and 2.8 eV above the VBM and separated by a 
Hund exchange splitting = 2 eV. The Cr e level con- 
tribution to the total magnetization is roughly 80 % of 
the total magnetization and the bonding ts states are 
weakly spin-polarized. The antibonding Iab levels are 
totally empty and forms two bound states around 4.3 eV 
above the VBM. The number of electrons located on the 
Cr ion is 22.3, which correspond to a Cr^+ configuration 
4s^3(i^ and spin Si = 1. Cr^+ induces two extra electrons 
distributed over NN carbon dangling bonds of configura- 
tion 2s'^2p^ and spin ^2 = 0. The Cr (^2)^ level hybridizes 
with the NN carbon 2j9^ t2 defect level. The bonding 
states are fully occupied and the e levels are partially oc- 
cupied in the configuration (e^)^(e^)^ corresponding to a 



FIG. 9. (color online) Top : Total, Ni and NN C partial DOS. 
Bottom : Schematic representation of the p-d hybridization 
model for Ni^. Symbols are the same as in Fig. [3| 

Fig. [9] shows the total, Ni and NN carbon partial 
SGGA density of states (DOS) of the NiCes supercell. 
The Ni 3d levels are split by the crystal field into e 
and t2 levels. Ni t2 levels hybridize with the NN car- 
bon 2p levels to form bonding ts and antibonding Iab 
hybrid levels. The spin up and spin down bonding Ib lev- 
els forms a crystal field resonance state situated around 
-2.1 eV below the VBM. The spin up and spin down 
e levels are localized on the Ni site, giving rise to two 
peaks approximately 1 eV below the valence band maxi- 
mum (VBM) and separated by a Hund exchange splitting 
= 0.4 eV. The spin up and spin down antibonding 
tAB levels are located at 2.5 and 3.0 eV above the VBM. 
The antibonding Iab level contribution to the total mag- 
netization is more than 70 %, with 50 % from Ni 3d level 
and 20 % from the NN carbon 2p level. The number of 
electrons located on the Ni ion is 26.9, which corresponds 
to a Ni+ configuration 4s^3d^ with a spin Si = ^. Ni+ 



8 



induces one extra electron distributed over the four NN 
carbon dangling bonds of configuration 2s^2p^ and spin 
S2 = \' The Ni ^2 levels hybridize with the NN carbon 
2p^ t2 defect level. This p-d hybridization picture iden- 
tifies the origin of the calculated SGGA DOS as due to a 
FM interaction between and 5*2 as shown on Fig. [9j 



The unusual case of Cu., 



I 3. S — ^ center. 



30 
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FIG. 10. (color online) Top : Total, Cu and NN C par- 
tial DOS. Bottom : Schematic representation of the p-d hy- 
bridization model for Cu^. Symbols are the same as in Fig.js] 



els are situated at 1.9 and 2.7 eV above the VBM. The 
number of electrons located on the Cu ion is 28.1, which 
corresponds to a Cu+ configuration As^'id^^ and a spin 
Si = 0. Cu+ induces one extra electron distributed over 
the NN carbon dangling bonds of configuration 2s^2p^ 
and spin 5*2 = |. The Cu ^2 level hybridizes with the 
NN carbon 2p^ t2 defect level. The 2p^ configuration of 
the NN carbon dangling bond is responsible for the total 
spin 5=1 {Mt = 3.0 /ie), which is distributed over 
the Cu (Mcu = 0.8 jub) and the 4 NN carbon atoms 
(Mc = 0.4 /xb)- Fig. [To] explains the observed DOS by a 



p-d hybridization model. Cu+ 3d levels are nearly spin 
degenerate corresponding to a non-magnetic 3d^^ closed- 
shell configuration. 



G. Outliers: Sc° (S = |) and Zn° (S = 1). 



Fig. 11 



shows the total, Sc and NN carbon partial 
SGGA density of states (DOS) of the ScCea supercell. Sc 
3d levels hybridize slightly with the NN carbon ^2 levels 
and give rise to 2 spin-up and 2 spin-down bound states 
situated at 0.7 and 1.2 eV (1.0 and 1.6 eV). The spin- 
down t2 level at 1.6 eV is partially empty (1 hole) which 
is responsible for the total spin-polarization Mt = 1 /xr 
which is distributed on the Sc and NN carbon atoms. The 
two spin-degenerated peaks at 4.6 eV arises from the Sc 
e levels. There are no antibonding levels and the density 
of states cannot be explained by the p — d hybridization 
model used above for Ti through Cu. 




Energy [eV] 



Fig. [To] shows the total, Cu and NN carbon partial 
DOS of the CuCes supercell. The Cu 3d levels are split 
by the crystal field into e and ^2 levels. Cu t2 levels 
hybridize with the NN carbon 2p levels to form bonding 
ts and antibonding Iab hybrid levels. The bonding ts 
levels are weakly spin-polarized and are located around 
2.6 eV below the VBM. The spin up and spin down e lev- 
els are localized on the Cu site, giving rise to two nearly 
spin-degenerate peaks situated around 2.4 eV below the 
VBM. The spin up and spin down antibonding Iab lev- 



FIG. 11. (color online) Total, Sc and NN C partial DOS. 

Fig. [12] shows the total, Zn and NN carbon partial 
SGGA density of states (DOS) of the ZnCea supercell. Zn 
gives rise to two bound states of t2 symmetry at 1.6 and 
2.1 eV above the VBM which are almost entirely of NN 
carbon 2p character. The partially empty (2 holes) bound 
state at 2.1 eV above the VBM is responsible for the total 
magnetization Mt = 2 ji^ which is distributed over the 
Zn and NN carbon atoms. The number of electron on Zn 
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is 29.0 corresponding to a configuration 3d^^4:S^. Thus 
the Zn 4s level is almost depleted and the Zn 45 electron 
occupies a bound state of mostly NN carbon 2p character 
to lower its energy. The Zn+ 3d levels are delocalized and 
are located at about -6 eV below the top of the valence 
band and are nearly spin degenerate, corresponding to a 
non-magnetic 3d^^ closed-shell configuration. 




FIG. 12. (color online) Total, Zn and NN C partial DOS. 



III. HIGH-SPIN TO LOW-SPIN TRANSITION 
WITH STRAIN 




fixed-spin calculations. The HS (LS) state of Cr^ and Ni^ 
corresponds to a total magnetic moment Mrs = 2 /iB 
(Mls = /ie) whereas the HS (LS) state of Cu^ corre- 
sponds to Mhs = 3 /iB (^LS = 1 Mb)- For Cr^, we do 
not obtain a transition due to the large magnetic energy 
of the 5* = 1 ground state, about 1 eV energetically lower 
than the S = non-magnetic solution. The energetic sta- 
bility of the magnetic ground state is associated with the 
e character of the highly spin-polarized states, whereas 
both Ni and Cu have magnetism driven by the t2 states. 
The results for Ni^ and Cu^ are presented in Fig. 13 
We obtain a transition from high-spin = 1 (5 = |) to 
low-spin S = {S = ^) under compressive hydrostatic 
strain with a transition at e// = — 7 % ( e// = — 6 %) for 
NiO (Cu^). 




-6 -4 
Hydrostatic Strain [%] 



FIG. 14. Low-spin state of Ni°, 5' = : Top : Total and Ni 
and NN C partial DOS. Bottom : Schematic representation 
of the p — d hybridization model for Ni^. 



FIG. 13. Total energy difference between low-spin and high- 
spin fixed spin moment SGGA calculations as a function of 
hydrostatic strain. 

We now study the effect of hydrostatic pressure on the 
magnetic properties of high-spin TM^ dopants (TM = 
Cr, Ni and Cu). For that purpose, we take the total en- 
ergy difference between high-spin (HS) and low- spin (LS) 



Fig. [T4]and[T5]presents the DOS of the low-spin SGGA 
ground state of Ni^ and Cu^ under compressive hydro- 
static pressure en = —10 %. The corresponding total 
magnetic moment for the (spin unfixed) SGGA solution 
is Mls = 0.1 /iB for Ni^ and Mls = 1-4 /iB for Cu^. The 
TM+ configuration 4s^3(i^^^ remains unchanged, which 
induces an extra electron on the NN carbon atoms, cre- 
ating a 2s^2p^ defect level. Fig. 14 and 15 explain the 
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observed DOS with a p — d hybridization model. For 
Ni^, the spin-degenerate Ni 3d levels hybridize with the 
spin-degenerate 2p^ dangling bonds giving rise to a non- 
magnetic interaction, corresponding to an antiferromag- 
netic interaction between the Ni /Si = ^ and a spin 
5*2 = distributed over the NN carbon dangling bonds. 
For Cu^, the spin-degenerate Cu 3d levels hybridize with 
the 2p^ dangling bond. The 2p^ configuration is respon- 
sible for the total spin S = ^ distributed over the NN 
carbon dangling bonds. 



— Total 

— Cu 




FIG. 15. Low-spin state of Cu^, 5* = | : Top : Total and Cu 
and NN C partial DOS. Bottom : Schematic representation 
of the p — d hybridization model for Cu^. 

The nearly-degenerate S = 1 and 5* = states of 
the substitutional nickel impurity in diamond can be un- 
derstood as an exchange-coupled system of two electron 
spins: one localized on the nickel ion and one delocal- 
ized on the four nearest-neighbor (NN) carbon atoms^^. 
Strain changes the exchange coupling and as a result 
modifies the ground state from parallel to antiparallel 
alignment. A similar phenomenon occurs for Cu, but the 
transition is between 5 = 3/2 and S = 1/2. The exis- 
tence of a pressure induced transition between a high-spin 



and low-spin state is also observed in transition-metal- 
ion compounds, where it was shown theoretically to be 
due to a competition between localization (induced by 
the Hund's exchange coupling which favors the high-spin 
state) and a tendency towards delocalization induced by 
the crystal field (which favors the low-spin state) 



IV. CONCLUSION 

We have studied by first-principle calculations the 
chemical trends of neutral substitutional transition-metal 
dopants in diamond. The origin of the magnetism for 
these dopants has been explained by a p-d hybridization 
model. Ti and Fe impurities are non-magnetic corre- 
sponding to closed bonding hybridized t2 and an empty 
(Ti) or full (Fe) e level respectively. For V, Cr and Mn, 
the magnetization is entirely localized on the TM ion and 
driven by the e levels. For Co, Ni and Cu, the magneti- 
zation is distributed over the TM ion and the NN carbon 
atoms with a magnetization due to the antibonding t2 
levels. Electron paramagnetic study of these magnetic 
centers would be of great interest in order to confirm the 
calculated chemical trend. 

The calculated total magnetic moments for the 
TM ion series are in agreement with previous GGA 
calculations^^. Ref. 'SF finds the same relaxation with 
point goup Td symmetry for Ti, Cr and Fe. They found 
a NN relaxation of D2d for V, Mn and Co and a Ci sym- 
metry relaxation for Ni. To verify the validity of our 
results, which were all- atom relaxation calculations per- 
formed with Td symmetry, we did nearest-neighbor atom 
relaxation calculations within GGA (1 meV/A precision) 
with no symmetry constraints, fixing the other atomic 
positions to their ideal value. For Ni, the tetrahedra cal- 
culated this way is slightly distorted to Ci symmetry, 
however the distortion is less than a 1% difference in dis- 
tance and less than a 0.1% difference in angle, with a 
magnetic energy Em = 48 meV and a total spin S = 1. 
Hence our T^^-symmetry results represent well the state 
of these transition-metal ions. 

The nearly-degenerate 5* = 1 and 5* = states of 
the substitutional nickel impurity in diamond can be un- 
derstood as an exchange-coupled system of two electron 
spins: one localized on the nickel ion and one delocalized 
on the four nearest-neighbor (NN) carbon atoms. Such 
nearly-degenerate S = 1 and S = states can be used 
to construct an effective two-state decoherence free sub- 
space from the 5'^ = triplet (Tq) and the singlet (S) 
states. This approach would be similar to that imple- 
mented for double quantum dots^^ with electrostatic gat- 
ing. For the nickel ion, strain modulation could be used 
to manipulate the energy splitting between the S and the 
To states, instead of the electrostatic gating used for dou- 
ble quantum dots^^. This is one way that the electronic 
configurations of transition-metal dopants could be used 
to encode quantum information in a form amenable for 
solid state quantum information processing. 
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